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(54) Procede de construction et de codage simple et systematique de codes Ldpc 



(57) Le procede permet la construction de codes 
LDPC comportant N symboles dont K libres chaque co- 
de etant defini par une matrice de contrdle A comportant 
M=IM-K lignes N colonnes et t symboles non nuls dans 
chaque colonne. II consiste: 

a attribuer a toutes les lignes de la matrice de con- 
trdle A un meme nombre de symboles non nuls, 
a prendre comme nombre de symboles t un nombre 



impair le plus petit possible, 
a definir les colonnes de facon que deux colonnes 
quelconques de la matrice de controle A n'aient au 
plus qu'une valeur commune non nulle, 
et a definir les lignes de fagon que deux lignes de 
la matrice de controle A n'aient qu'une valeur com- 
mune non nulle. 

Application: Transmissions radio. 
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Description 

[0001] La presente invention concerne un procede de construction et de codage simple et systematique de codes 
connus sous I'abreviation anglo-saxonne LDPC de Tow Density Parity Check". 
5 [0002] Les codes de Gallager, proposes vers 1 963, sont a rorigine des codes LDPC actuellement envisages comme 
alternative aux turbo-codes. 

[0003] Un article publie dans la revue IEEE Transaction on Information Theory Vol 45 n°2 Mars 1999 de MJ C 
MacKay ayant pour titre "Good Error Correcting Codes Based on very Sparse Matrices" presente des resultats inte- 
ressants concernant ces codes, en particulier le fait que: 

10 

ce sont asymptotiquement pour des blocs de taille elevee de "tres bons codes" 

le decodage pondere ("soft decoding", ou "decodage souple") est facile a mettre en oeuvre. 

[0004] Mais il n'existe pas de methode autre qu'heuristique pour les construire. 
15 [0005] Suivant cette technique de codage un code (N, K) comportant N symboles dont K libres est defini par sa 
matrice de controle de parite A, comportant M = N - K lignes et N colonnes. 

[0006] La matrice de controle A est caracterisee par sa faible "densite" : il faut entendre par la qu'elle comporte un 
nombre reduit d'elements non nuls. 

[0007] Plus precisement, elle comporte exactement t symboles non nuls dans chaque colonne, tous les autres etant 
20 egaux a 0. 

[0008] Si les symboles d'un mot de code sont notes ci, i = 0...N-1 et les elements de la matrice de controle Aij, le 
code satisfait M = N-K relations de la forme : 

25 z /-o...N-i A ml c iP ourm=0...M-1 

[0009] Les methodes proposees par M J C MacKay consistent a batir une matrice A initiale a partir de matrices unite 
ou tridiagonales plus petites, puis a permuter leurs colonnes pour arriver au resultat souhaite. L'experience montre 
cependant qu'il est difficile de satisfaire les differentes contraintes imposees pour leur construction. 
30 [0010] Le but de invention est de pallier les inconvenients precites. 

[001 1] A cet effet invention a pour objet, un procede de construction de codes LDPC comportant N symboles dont 
K libres chaque code etant defini par une matrice de controle A comportant M=N-K lignes N colonnes et t symboles 
non nuls dans chaque colonne caracterise en ce qu'il consiste : 

35 a - a attribuer a toutes les lignes de la matrice de controle A un meme nombre de symboles T non nuls, 

b - a prendre comme nombre de symboles T un nombre impair le plus petit possible, 

c - a definir les colonnes de facon que deux colonnes quelconques de la matrice de controle A n'aient au plus 
qu'une valeur non nulle, 

d - et a definir les lignes de facon que deux lignes de la matrice de controle A n'aient qu'une valeur commune non 
to nulle. 

[0012] Le procede selon I'invention a pour avantage qu'il permet de simplifier les algorithmes de codage et de de- 
codage, en utilisant une matrice de controle A qui est la moins dense possible tout en donnant de bonnes performances 
pour une complexity raisonnable, la puissance de calcul necessaire etant proportionnelle au nombre t. Dans la mesure 
45 ou il y a peu d'erreurs, la contrainte "c" ci dessus permet a I'algorithme de decodage de toujours converger. 

[0013] D'autres caracteristiques et avantages de I'invention apparaitront a I'aide de la description qui suit faite en 
regard des dessins annexes qui represented: 

[0014] La figure 1 un tableau representant une partition de la matrice de controle A. 

[0015] Les figures 2 et 3 des groupements des m 2 sous matrices de gauche du tableau de la figure 1 en une sous 
so matrice MxM et n-m sous-matrices MxP 

[0016] Les figures 4 et 5 une matrice de controle A obtenue respectivement selon une premiere et une deuxieme 
variante d'execution du procede selon I'invention. 

[0017] Les figures 6 a 9 une matrice de passage pour le calcul des symboles de redondance. 
[001 8] Pour la mise en oeuvre du procede selon I'invention la matrice de controle, A, est subdivisee comme le montre 
55 la figure 1 en n sous-matrices de M lignes et P colonnes ou m x n sous-matrices carrees de P lignes et P colonnes 
telles que N=nP et M=mP et n et m premiers entre eux. 

[001 9] Les m 2 sous-matrices de gauche sont ensuite regroupees comme le montre la figure 2 en une sous-matrice 
MxM (ceci permettra de simplifier grandement ralgorithme de codage) et les autres en n - m sous-matrices MxP. 
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[0020] Le procede de construction est decrit ci - apres selon deux variantes, selon que m est egal a 1 ou t. 
[0021] Les valeurs differentes de m ne sont pas a envisager ici, a cause de ta condition "a" qui exige que tr soit 
entier. En effet, tr = t N / M, soit encore t n / m. 

[0022] n et m etant premiers entre eux, m doit diviser t et ne peut done etre qu'egal a 1 ou t pour t premier et petit 

5 (vrai pour les faibles valeurs de t, a savoir 3, 5, 7). 

[0023] Dans la premiere variante ou m = 1 (codes de redondance r/ (r - 1)) f le procede selon I'invention est valable 
pour les codes de taux de redondance N/K de la forme r / (r - 1), ou le nombre de symboles de redondance est 
exactement N / r (r entier). Dans ce cas, M est egal a P et le tableau de la figure 2 se resume au tableau de la figure 
3. Le procede consiste alors a chercher n sequences de longueur M comportant t "1" et (M - 1) "0". 

10 [0024] Ces sequences notees ci-apres w[0..n-1], sont obtenues par : 

une autocorrelation cyclique egale a 0, 1, ou t (la sequence i decalee ne coincide avec elle-meme non decalee 
qu'en 0 ou 1 point) telle que pour tout i = 0...n-1 : 

15 

£ /c=o....m-i wMM wflPl = 1 (P ar definition) 



20 



25 



z k=o..m-i w Wt k ] w [']K k + P) modulo M] = 0 ou 1, pour p = 1..M-1 

une intercorrelation cyclique egale a 0 ou 1 (la sequence I decalee ou non ne coincide avec la sequence j qu'en 
0 ou 1 point) telle que pour toute paire {i, j} ou i = 0...n-1 et j = 0...n-1 sont differents : 

z k=o....am w I'JM w DlK k + P) modulo M] = 0 ou 1, pour p = 0..M-1 



[0025] L'algorithme de calcul des sequences w est tres simple : ii determine successivement les positions pos[0][0... 

t-1], pos[1][0...t-1] pos[n-1][0...t-1], ou ces sequences possedent un "1", en commen* ant par pos[x][0] = 0, pos[x] 

[1] = 1, pos[x][t-1] = t-1, et en les modifiant pour satisfaire les conditions d'autocorrelation et intercorrelation. 
30 [0026] Pour t = 3, les boucles de calcul mises en oeuvre sont montrees a I'annexe 1. 

[0027] Cet algorithme echoue lorsque n est trop grand compte tenu de M : peu de "petits" codes sont trouves, mais 
ceci est de peu d'importance car on cherche generalement des codes de grande taille (N » 100). 
[0028] Les colonnes de la matrice A sont alors tout simplement les vecteurs w permutes circulairement : 

35 - kerne sous-matrice (k = 0...n-1) 

A[ligne][colonne] = w[k][(ligne - (colonne - k P)) modulo M] 
Avec: 

ligne = 0...M-1 

et colonne = kP...(k + 1)P - 1 

40 

[0029] Ainsi, chaque ligne de A comporte exactement t valeurs non nulles dans chacune des n sous-matrices, soit 
un total de n = tr. 

[0030] Un exemple de matrice A, obtenue par ce procede pour un code LDPC(75, 50, t=3, tr=9) de redondance 3 / 
2 (r = n = 3) avec P = 25 est montre a la figure 4. Selon le tableau represents on peut constater que: 

45 

w[0][i] = 1 pour i = 0, 1, 3 
w[1][i] = 1 pouri = 0, 4,9 
w[2][i] = 1 pour = 0, 6, 13 

so [0031] La construction proposee garantit que : 

chaque colonne comporte exactement t valeurs non nulles (par definition des w) 
- chaque ligne comporte exactement tr valeurs non nulles (grace aux proprietes d'autocorrelation et intercorrelation 
des w) 

55 . toute paire de colonnes distinctes a au maximum une valeur non nulle commune (idem) 
toute paire de lignes distinctes a au maximum une valeur non nulle commune (idem) 

[0032] Selon une deuxieme variante inspiree de la precedente correspondent au cas ou m=t le procede selon I'in- 
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vention recherche n - m + 1 sequences de longueur M comportant t "1" et (M - 1) "0", sequences notees w[0..n-m]. 
[0033] La premiere sequence, w[0] est obtenue par : 

- une autocorrelation cyclique egale a 0, 1 , ou t (la sequence 0 decalee ne coincide avec elle-meme non decalee 
qu'en 0 ou 1 point) telle que: 

*K=o....M-iW[0][k] w[0][k] = t (par definition) 

£k=o....m-i w [°]M w0][(k + p) modulo M] = 0 ou 1, pour p = 1..M-1 

[0034] En fait, c'est la meme definition que pour m = 1 . 
[0035] Les sequences suivantes, w[1 ...n-m] sont obtenues par 

une autocorrelation cyclique egale a 0 ou t pour des decalages multiples de m (la sequence i decalee d'un multiple 
de m ne coincide jamais avec elle-meme non decalee) telle que: 
pour tout i = 1...n-m : 

z k=o....am w H[k] w[i][k] = t (par definition) 

£ k=o...m-i w [i][k] w[i][(k + p m) modulo M] = 0, pour p = 1..P-1 

par une intercorrelation cyclique egale a 0 ou 1 avec la sequence w[0] (la sequence i decalee ou non ne coincide 
avec la sequence 0 qu'en 0 ou 1 point) telle que: 
pour tout i = 1...n-m : 

£ k=o...w-i w ['H k l w I°][( k + P) modulo M] = 0 ou 1, pour p = 0..M-1 



- et par une intercorrelation cyclique avec les sequences w[1 ..n-m] egale a 0 ou 1 pour des decalages multiples de 
m (la sequence i decalee ou non d'un multiple de m ne coincide avec la sequence j qu'en 0 ou 1 point) telle que: 
35 pour toute paire {i, j} ou i = 1 ...n-m et j = 1 ...n-m sont differents : 

z /c=o....m-i wfllM w D][(k + pm) modulo M] = 0 ou 1, pour p = 0..P-1 

40 [0036] L'algorithme de calcul des sequences w est le meme que precedemment. Seuls changent les criteres d'auto- 
correlation et d'intercorrelation, celles-ci n'etant a verifier que sur P points au lieu de M. 

[0037] Les colonnes de la matrice A sont alors les vecteurs w permutes circulairement avec un pas egal a 1 ou m 
telles que: 

45 - Sous-matrice MxMa gauche de A : 

A[ligne][colonne] = w[0][(ligne - colonne) modulo M] 
Avec: 

ligne = 0...M-1 
colonne = 0...M - 1 

50 

Sous-matrices MxP suivantes (en nombre egal a n - m) pour k = m...n-1 : 
A[ligne][colonne] = w[k - m + 1][(ligne - m (colonne - k P)) modulo M] 

Avec: 

ligne = 0...M-1 
55 colonne = kP...(k + 1)P - 1 

[0038] Ainsi, chaque ligne de A comporte exactement m = t valeurs non nulles dans ses M premieres colonnes, puis 
1 valeur non nulle dans chacun des n-m paquets de P colonnes successives, soit un total de n ou tr. 
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[0039] Un exemple de matrice A, obtenu selon la deuxieme variante du procede selon Tinvention est montre a la 
figure 5 pour un code LDPC(75, 30, t=3, tr=5) de redondance 5 / 2 (n = 5, m = 3) avec P = 15. On constate sur le 
tableau de la figure 5 que: 

5 w[0][i] = 1 pour = 0, 1,3 

w[1][i] = 1 pouri = 0, 4, 8 
w[2][i] = 1 pour = 0, 5, 10 

[0040] Le procede selon I'invention qui vient d'etre decrit sous ses deux variantes conduit directement a un algorithme 
10 de codage de symboles de redondance Yi et de symboles d'information Xi tres simple. 

[0041] Pour cela, il suffit de considerer que les symboles de redondance Yi sont les M premiers symboles d'un mot 
de code, et les symboles libres Xi (information) sont les N - M derniers. 

[0042] Les equations que doivent verifier tout mot de code peuvent done etre reecrites sous la forme : 

15 

£ /=o...m-i Ami Yi + s /=m.../v-i Ami Xi = °« P° ur m=0...M-1 

ou encore : 

20 

£ /=o...m-i Ami Yi = Zm , pour m=0...M-1 

avec 

25 

Zm = - S /!=M w-1 Ami Xi, pour m=0...M-1 

[0043] Le procede consiste alors a calculer dans un premier temps les M quantites Zm de la matrice de passage, 
et ensuite les symboles de redondance : 

30 

Ym = £, =0 ...AM Bmi Zi » P° ur m=0...M-1 

[0044] A titre d'exemple, pour le code LDPC(75, 50) les quantites Zm sont calcules par le systeme d'equation defini 
35 par le tableau de la figure 6 qui, apres resolution, se transforme en le tableau des symboles de redondance de la figure 7. 
[0045] La matrice B d'element generique Bij est ('inverse de la partie gauche (de dimension MxM)de la matrice 
A. Elle a une forme tres simple : par construction, toutes ses colonnes sont des permutations circulaires de la sequence 
w[0][0..M-1] : 

Aij = w[0][(i - j) modulo M], i=0...M-1, j=0...M-1 
40 [0046] La matrice B comporte alors M lignes qui sont des permutations circulaires d'une ligne unique b[0...M-1], a 
savoir : 

Bij = b[(j - i) modulo M] 

B etant I'inverse de P^, les coefficients b sont definis par : 
45 x /=o.../w-i w[0][i] b[(i + k) modulo M] = 1 si k = 0, 0 si k = 1 ...M-1 

[0047] Par exemple, pour le code LDPC(75, 50), les coefficients de redondance Ym sont calcules par le systeme 

d'equations defini par le tableau de la figure 8 qui apres resolution se transforme en le tableau de la figure 9. 
so [0048] Cependant, il existe des cas ou le calcul est impossible. 

[0049] On peut en effet ecrire les equations qui les definissent sous la forme : 
,a m t{b[0], b[1], .... b[M-1]} = t{1, 0, 0, 0} 

[0050] La matrice *A M est une matrice circulante, dont la premiere ligne est egale a a[0...M-1] = w[0]. 

[0051] Son determinant est egal au produit de ses M valeurs propres I0...M-1 . 
55 [0052] La kerne valeur propre est elle-meme donnee par : 

= S /=0...M-1 a Pl a ik 
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ou a est une racine Meme de I'unite. 
[0053] Par exemple: 

- pourw[0] = a[0...M-1] = {1, 1,0, 1,0,0,0,....} 

- pour des codes binaires (on se situe dans le corps de Galois CG(2) ou I'addition est equivalente au OU EXCLUSIF 
(XOR) et la multiplication au ET logique) 

on a : 

X 1 =1 + a + a 3 
[0054] Si M est multiple de 7, il se trouve que Tequation 



definit un corps de Galois ou a est racine 7eme de I'unite (le polynome g(x) = 1 + x + x 3 est irreductible et primitif 
dans CG(2) et genere un corps de Galois CG(23)), ce qui signifie que ^ = 0. 
20 [0055] Parmi les codes LDPC trouves par ralgorithme propose, il faut done eliminer ceux ou M est multiple de 7 si 
Ton garde ce w[0] ci car : 

Tune des valeurs propres de i A fA sera nulle 
done son determinant sera nul 
25 - done on ne pourra pas trouver de b[i] convenables 

done on ne pourra pas effectuer le codage (calculer les Yi) 

[0056] D'une maniere tres generate, quel que soit le choix fait pour w[0], il y aura des valeurs de M ne convenant 
pas car ne permettant pas de faire le codage. 
30 [0057] On montre facilement (en factorisant x M - 1 et a(x) = Z f -JlM a[i] x*) que ces valeurs de M sont les multiples 
d'une valeur M0 pour laquelle a(x) divise xMO - 1 . 
[0058] Par exemple, pour des codes binaires avec t = 3: 

- w[0] = {1,1,0,1,...} 
35 - w[0]={1,0, 1,1,...} 

interdisent M multiple de 7 (a(x) definit une racine 7eme de I'unite) 

- w[0] ={1,1, 0,0,1,...} 
40 - w[0] = (1,0,0, 1,1,...} 

interdisent M multiple de 15 (a(x) definit une racine 15eme de I'unite) 

- w[0] = {1, 0, 1, 0, 1, ...} n'est pas accepte (autocorrelation incorrecte) 
45 - w[0]= 0.1,0,0,0,1,...} 

- w[0]={1.0,0,0,1,1,...} 

interdisent M multiple de 3 (a(x) est multiple de 1 + x + x2 qui definit une racine 3eme de I'unite) 

so - w[0]= {1,0, 1,0,0,1,...} 

- w[0] ={1,0, 0,1,0, 1,...} 

interdisent M multiple de 31 (a(x) definit une racine 31eme de I'unite) 
[0059] Le calcu! des coefficients b[i] est effectue de la maniere suivante: 
55 [0060] Pour une valeur de M non interdite, il existe un algorithme particulierement simple de calcul des b[i] a partir 
des a[i] (ou w(0][0...M-1]). Cet algorithme repose sur I'observation que la serie des b[i], apres periodisation et filtrage 

par un filtre a reponse impulsionnelle finie (RIF) A(z) dont les coefficients sont les a[M-1, M-2 1, 0] doit donner la 

serie {1, 0, 0, ...} periodisee. En fait, pour un code binaire utilisant I'un des w[0] enumeres precedemment, cette serie 
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est formee de la concatenation de sequences maximales (Maximal Length Sequences) de longueur 7 (ou 15 ou 31 
ou 63). 

[0061] On calcule done la reponse impulsionnelle du filtre a reponse impulsionnelle infinie (Rll) 1 / A(z) et on en 
extrait une tranche de longueur M qui, une fois periodisee, donne la serie {1 , 0, 0, ...} apres filtrage par A(z). 
5 [0062] Par exemple, pour un code binaire avec t = 3 et pour lequel seuls a[0], a[k1] et a[k2] ne sont pas nuls, un 
algorithme correspondant est fourni a I'annexe 2. 

[0063] Dans un esprit de simplification pour ne pas effectuer le calcul precedent a chaque codage, Talgorithme de 
codage peut encore etre defini par les k2 derniers elements b[M-K2...M-1] qui, par recurrence (filtrage par 1 / A(z)) 
permettent de recalculer tous les autres. 

10 [0064] Egalement comme la deuxieme phase de Palgorithme de codage standard (calcul des Y a partir des 2) com- 
porte en moyenne M 2 / 2 operations, ce qui peut devenir important pour des codes consequents : la complexity etant 
une fonction quadratique de la taille, et que de plus, il est necessaire de stocker le tableau intermediate Z (M elements) 
et connaitre le tableau b (M elements aussi) s'il n'est pas calcule sur place, cette partie de l'algorithme peut etre modifiee 
pour n'utiliser que deux tableaux intermediates de tres petite taille en reecrivant ainsi les equations donnant les Y de 

15 la facon montree par le tableau de la figure 9 (exemple pour un code LDPC(75, 50)) : 

[0065] Les M - k2 (pour t = 3) premieres lignes sont les M - k2 dernieres lignes du systeme d'equations donnant Y, 
avant resolution. 

[0066] Les k2 dernieres lignes sont les k2 dernieres lignes du systeme d'equations donnant Y, apres resolution. 

[0067] II suffit alors de calculer les Y dans I'ordre inverse, a savoir Y[M-1], Y[M-2] Y[0]. 

20 [0068] Le nombre d'operations a effectuer est alors en moyenne de k2 M / 2 (calcul de Y[M-1] ... Y[M-k2]) suivi de 
t(M - k2) (calcul de tous les autres) soit approximativement (t + k2/2) M : la complexity n'est plus qu'une fonction lineaire 
de la taille. 

[0069] L'algorithme utilise X[M...N] en entree. 

[0070] La partie basse de X (X[0...M-1]) est utilisee comme stockage temporaire pour les Z : X[0...M-1] stocke Z 
25 [k2...M-1 , 0...k2-1] pour eviter un decalage circulaire en phase finale. 

[0071] Les b[i] sont calcules iterativement sur place, a partir des b[M-k2...M-1]. 
[0072] Le code est defini par deux tableaux : 

- le tableau finB[0...k2-1] des k2 derniers elements de b 

30 - le tableau pos[0...(n - m + 1)t] contenant les positions des elements non nuls des sequences w[0], w[1], .... w[n-m]. 

On utilise deux buffers internes de taille k2 : 

- reg[0...k2-1] pour calculer les b[i] 

35 - temp[0...k2-1] pour stocker les valeurs intermediates de Y[M-k2...M-1]. 

[0073] L'algorithme complet de codage rapide est alors celui montre a Tannexe 3. 

[0074] Ces algorithmes sont tres simples a mettre en oeuvre. lis ont entre autre la caracteristique de definir un code 
par tres peu de parametres, a savoir les (n - m + 1)(t - 1) positions non nulles des "1" dans les sequences w et 
40 eventuellement k2 coefficients de codage. Meme s'ils ne donnent pas tous les codes possibles repondant aux condi- 
tions a-d (par exemple, pas le code (150, 125) de redondance 6/5, qui necessite n=6 sequences w de longueur P=25), 
its en donnent suffisamment pour que, dans une application quelconque ou N et K sont definis a priori, on puisse 
trouver soit 

45 - un code (NLDPC, KLDPC) avec NLDPC = N et KLDPC = K 

- un code voisin (NLDPC+d, KLDPC+d), avec d faible, qui sera raccourci par non transmission de d symboles utiles 
mis arbitrairement a zero. 

[0075] Par exemple, pourobtenirun codede (N, K) redondance 5/3 (taux 0.6), il suffit de partir d'un code (NLDPC+d, 
50 KLDPC+d) de redondance 8 / 5 (taux 0.625) avec d = NLDPC / 15. Pour des valeurs de N inferieures ou egales a 500 
et t = 3, il est possible de construire tres rapidement 932 codes differents dont les redondances sont les suivantes (on 
s'est volontairement limite aux redondances comprises entre 4 et 8/7, et aux codes ou w[0] = {1, 1, 0, 1, 0, 0, 0, ....} 
pour lesquels k2 = 3) : 



R = 


4/1 


soit 4.000 


(105 codes) 


R = 


5/2 


soit 2.500 


(82 codes) 


R = 


6/3 ou 2/1 


soit 2.000 


(203 codes) 
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(suite) 
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R 


= 7/4 


soit 1.750 


(55 codes) 


R 


= 8/5 


soit 1.600 


(47 codes) 


R 


= 9/6 ou 3/2 


soit 1.500 


(124 codes) 


R 


= 10/7 


soit 1.428 


(34 codes) 


R 


= 11/8 


soit 1.375 


(28 codes) 


R 


= 12/9 ou 4/3 


soit 1.333 


(84 codes) 


R 


= 13/10 


soit 1300 


(20 codes) 


R 


= 14/11 


soit 1.273 


(17 codes) 


D 

r\ 


- lonz ou o/4 


SOU \.ZD\J 


(Ob cooes) 


R 


= 16/13 


soit 1.231 


(11 codes) 


R 


= 17/14 


soit 1.214 


(7 codes) 


R 


= 18/15 ou 6/5 


soit 1.200 


(34 codes) 


R 


= 19/16 


soit 1.187 


(3 codes) 


R 


= 20/17 


soit 1.176 


(2 codes) 


R 


= 21/18 ou 7/6 


soit 1.167 


(17 codes) 


R 


= 24/21 ou 8/7 


soit 1.143 


(3 codes) 



[0076] De plus, toujours pour une valeur donnee de N inferieure ou egale a 500, il peut y avoir jusqu'a 12 codes 
differents (pour N = 480). 

[0077] Par exemple, des que N est multiple de 6 et superieur ou egal a 288, il existe toujours trois codes de longueur 
25 N et de redondances 6/5, 3/2 et 2/1 , par exemple LDPC(300, 250) + LDPC(300, 200) + LDPC(300, 150). 

[0078] Ceci est tres utile pour proteger efficacement un train binaire forme de trois trains binaires chacun de longueur 
N et de sensibilites differentes. 

[0079] Bien entendu, il est toujours possible d'envisager de nombreuses variantes de ces algorithmes, comme par 
exemple une permutation aleatoire des lignes et/ou colonnes de la matrice A. 
30 [0080] II est aussi important de signaler que I'adaptation a des codes non binaires est particulierement simple. 

ANNEXE 1 

35 for(x=0; x<n - m + 1 ; x++) { 

pos[x][0] = 0; 

for(pos[x][1]=pos[x][0]+1; pos[x][1]<M-1; pos[x][1]++) { 

for(pos[x][2]=pos[x][1 J+1 ; pos[x][2]<M; pos[x][2]++) { 
(si les conditions ne sont pas satisfaites, 
continuer 
sinon, aller a ok) 

} 

i 

(arreter : impossible de trouver un choix convenable 
so pourpos[x][0...t-1]) 

ok:; 

} 

55 
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ANNEXE 2 

(langage C : I'operateur " A " correspond au OU EXCLUSIF. 
/* Initialisation du passe de b, de longueur M 7 
for(i=M-k2; i<M; i++) 
b[i] = 0; 

/* Calcul de N valeurs successives de la 
reponse impulsionnelle de 1/A(z) 7 
b[0] = 1; 

for(i=1; kk2; i++) 

b[i] = b[(i+M-(k2-k1)) % M] A b[(i+M-k2) % M]; 
for(i=k2; i<M; i++) 

b[i] = b[i-(k2-k1)]Ab[i-k2]; 

I* S'arranger pour qu'il n'y ait qu'un 1 

dans les k2 dernieres positions de b filtre par A(z) 7 

weight = 0; /* tout sauf 1 7 

while(weight != 1) { 

/* Decaler d'un cran 7 
for(i=1; i<M; i++) 

b[M] = b[i]; 
b[M-1] = b[M-1-(k2-k1)] A b[M-1-k2]; 
/* Verifier 7 
weight = 0; 
for(i=M-k2; i<M; i++) { 

char sum = b[i] a b[(i+k1 ) % M] a b[(i+k2) % M]; 

if(sum) { 

shift = M - i; 
weight++; 

} 

} 
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r Cas particulier ou M est interdit 7 
if (weight == 0) 

retum(M_FORBIDDEN); 

} 

/* decalage circulaire final a droite : 
b[i] = b[(i - shift) % M] 7 
for(dec=0; dec < shift; dec++) { 

char temp = b[M-1]; 

for(i=M-1;i>0; i-) 
b[i] = b[i-1]; 

b[0] = temp; 

} 

retum(OK); 
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ANNEXE 3 

(langage C) : 

/* Phase 1 : calcul des M parites intermediates z. 
Ces parites sont calculees en lisant les colonnes successives de la 
matrice de codage, a savoir A[*][M], .... A[*][N] 

Elles sont mises en tete de x a titre temporaire 7 

#define z x 
for(i=0; i<M; i++) 
z[i] = 0; 

/* Boucle sur les n-m sous-matrices a la droite de A 7 
cO = M; 
d = cO + P; 

for(k = 1 ; k <= n - m; k++) { 
offset = 0; 

for(c = cO; c < d ; C++) { 
if(x[c] !=0) 

for(i=0; i<t; i++) { 

r p devrait etre offset + pos[i]. 

On le decremente de k2 pour eviter 

le decalage 

du tableau z avant la phase 3 V 
p = offset + pos[k*t + i] - k2; 
if(p<0) 

z[p + M] = z[p + M] M; 

else 

if(p < M) 

2[P] = Z[ P ]A1; 

else 

z[p - M] = z[p - M] M; 
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} 

offset = offset + m; 

} 

cO = d ; 
d =c1 + P; 

} 

/* Phase 2 : calcul des k2 derniers symboles de parite 7 
ixbO = M - 1 - k2; 

/* 1 : initialisation des k2 derniers elements de y 
temp[0...k2-1] = y[M-1, M-2, ...M-k2] 7 
for(k=0; k<k2; k++) 
temp[k] = 0; 

/* 2 : recopie des k2 derniers elements de b 
reg[0...k2-1] = b[M-k2...M-1] 7 
for(i=0; i<k2; i++) 

reg[i] = finB[i]; 

/* 3 : calcul iteratif des k2 derniers symboles 7 
for(i=0; i < M; i++) { 

r b[i] = {100 ...} * b[i-(k2-k1)] * b[i-k2] 

avec b[i-k2]...b[i-1] = reg[0...t2-1] 

On doit verifier : 

b[-k2] + b[k1-k2] + b[0] = 0 

b[-2] + b[k1-2] + b[k2-2] = 0 
b[-1] + b[k1-1] + b[k2-1] = 0 
b[0] + b[k1] + b[k2] = 1 
b[1] + b[1+k1] + b[1+k2] = 0 

...7 

if (i == k2) 

input = 1; 
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else 

input = 0; 
bi = input A reg[0] A reg[k1]; 
for(k=1;k<k2; k++) 

reg[k- 1] = reg[k]; 
reg[k2-1] = bi; 



If (bi != 0) 

« for(k=0; k<k2; k++) 



20 



25 



30 



35 



40 



if(z[(ixb0 - k + M) % M] != 0) 
temp[k] = tempfk] A 1 ; 

ixbO = ixbO + 1 ; 
if(ixbO == M) 
ixbO = 0; 



} 



I* 4 : Les z ont deja ete decales a gauche pour 
eviter Tecrasement. Sinon, il faudrait faire : 

for(k=0; k<M - k2; k++) 
z[k] = z[k + k2]; 

Recopie de temp a la fin de y 7 
#define y x 
for(k=0; k<k2; k++) 

y[M - 1 - k] = temp[k]; 



r Phase 3 : calcul de y[M-k2-1, M-k2-2, 0] 
*s y[k + k2 - k2] + y[k + k2 - k1] + y[k + k2 - 0] + z[k + k2] = 0 

y[k] va en x[k] 
z[k + k2] est en x[k] 

50 DOnC : 

x[k + k2 - k2] + x[k + k2 - k1] + x[k + k2 - 0] + x[k] = 0 
Soit : 

x[k + k2 - k2] = -(x[k + k2 - k1] + x[k + k2] + x[k]) 

55 
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10 



7 

for(k = M-k2-1;k>=0; k«) 

y[k] = y[k + k2 - k1] A y[k + k2] A z[k]; 



Revendications 



1. Procede de construction de codes LDPC pour la protection de trains binaires d'information, chaque train etant 
compose de N symboles decomposes en N-M symboles d'informations utiles Xj et M symboles d'informations de 
*5 redondance Y m , chaque code etant defini par une matrice de controle A comportant M=N-K lignes de N colonnes 

et t symboles non nuls dans chaque colonne, caracterise en ce qu'il consiste : 

a attribuer a toutes les lignes de la matrice de controle A un meme nombre de symboles non nuls, 
a prendre comme nombre de symboles t un nombre impair le plus petit possible, 
20 a definir les colonnes de facon que deux colonnes quelconques de la matrice de controle A n'aient au plus 

qu'une valeur non nulle, 

et a definir les lignes de facon que deux lignes de la matrice de controle A n'aient qu'une valeur commune 
non nulle. 

25 2. Procede selon la revendication 1 caracterise en ce qu'il consiste: 

a subdiviser la matrice de controle A de M lignes et N colonnes en n sous matrice de M lignes et P colonnes 
pour former n sous matrices de P lignes et P colonnes et a regrouper les m 2 sous matrices de gauche en une 
sous matrice MxM et les autres en n-m sous matrice MxP, 
30 - et a determiner Mxn sequences de vecteurs colonne w[0....n-1] de longueur M comportant t valeurs non nulles 

et (M-t) valeurs nulles en effectuant une autocorrelation et une intercorrelation cyclique des vecteurs colonne w. 

3. Procede selon la revendication 1 caracterise en ce qu'il consiste : 

35 - & subdiviser la matrice de controle A de M lignes et N colonnes en n sous matrices de M lignes et P colonnes 

pour former nxm sous matrices de P lignes et P colonnes et a regrouper les m 2 sous matrices de gauche en 
une sous matrice MxM et les autres en n-m sous matrice MxP, 

a determiner n-m+1 sequences de vecteurs colonnes w[0....n-m] de longueur M comportant t valeurs non 
*o nulles et (M-t) valeurs nulles, 

la premiere sequence w[0] etant obtenue par autocorrelation cyclique egale a 0, 1, ou a la valeur t de 
facon que la sequence w[0] decalee ne coincide avec elle meme non decalee qu'en 0 ou 1 point, 
les n-m sequences suivantes w[i][k] etant obtenues : 

45 - par une autocorrelation cyclique de valeur nulle ou egale a la valeur t de facon que la valeur de la 

sequence w[i] decalee d'un multiple de m ne coincide jamais avec elle-meme non decalee.et par une 
intercorrelation cyclique 

et par une intercorrelation cyclique de valeur nulle ou 1 avec les sequences w[1....n-m] pour des 
decalages multiples de m de facon que une sequence i decalee ou non d'un multiple de m ne coincide 
50 avec la sequence j qu'en 0 ou 1 point. 

4. Procede selon Tune quelconque des revendications 1 a 4 caracterise en ce qu'il consiste a determiner pour le 
codage d'informations utiles Xj une matrice de passage Z m egale au produit de la matrice de controle A par un 
vecteur colonne representant N-M symboles d'information Xj et a adjoindre aux symboles d'information des sym- 

55 boles de redondance Y m obtenus en resultat du produit de la matrice de passage Z m par une matrice B egale a 

Tinverse de la partie de dimension MxM de la matrice de controle A. 
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